Electric dipole rovibrational transitions in HD molecule 
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The rovibrational electric dipole transitions in the ground electronic state of the HD molecule are studied. A 
simple, yet rigorous formula is derived for the transition rates in terms of the electric dipole moment function 
D(R), which is calculated in a wide range of R. Our numerical results for transition rates are in moderate 
agreement with experiments and previous calculations, but are at least an order of magnitude more accurate. 

PACS numbers: 33.70.Ca, 95.30.Ky, 31.15.-p, 31.30. -i 



I. INTRODUCTION 

The electric dipole rovibrational transitions in the HD 
molecule are possible due to different masses of the proton 
and of the deuteron and thus slightly different binding ener- 
gies in hydrogen and deuterium. These transitions, for the 
first time, were observed by HerzbergjhjLand since then mea- 
sured by several groups 

ElSS&LZlIi- Theoretical 
calculations of the dipole transition moment were first car- 
ried out by Wick ifToll . somewhat later by Wu ifTHl and by 
Blinder yj . More elaborate calculations include those of 
Bunker lll3ll . Wolniewicz J3, Ford and Brow ne lfl5ll : and 
Thorson et al. d]. The most recent works ifll Jl5j [l6|l 
are in generally good agreement with experimental results in 

In this work we derive a compact formula for the dipole 
transition moment using a unitary transformation of the 
Hamiltonian followed by the adiabatic approximation, and 
present results in terms of the electric dipole moment func- 
tion D{R). We obtain D(R) for a wide range of internu- 
clear distances R 6 (0.5, 12) au, which enables calculations 
of electric dipole transitions between all rovibrational states 
including the highly excited ones. Although they have not 
been measured, these dipole transitions between highly ex- 
cited states together with electric quadrupole transitions lead 
to the cooling of the hydrogen clouds in the interstellar space 
lfl7tl . The obtained transition rates between low lying rovi- 
brational states are the most accurate to date, and agree with 
experimental values with minor exceptions. 



II. DERIVATION OF THE TRANSITION DIPOLE 
MOMENT 

In order to derive a formula for the dipole transition mo- 
ment, we consider a diatomic molecule in the reference frame 
of the geometrical center of the two nuclei. The total wave 
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function is a solution of the stationary Schrodinger equation 

H0 = E<p, (1) 

with the Hamiltonian 

H = H cl + H n , (2) 

split into the electronic and nuclear parts. In the electronic 
Hamiltonian 



E 



2 m c 



V. 



(3) 



with V including the Coulomb interaction, the nuclei have 
fixed positions Ra (proton) and Rb (deuteron), while the nu- 
clear Hamiltonian is 

\ r 1% \ ., ' ~ 



V 2 



(Eg Va ) 
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where R — Ra — Rb and /x n is the nuclear reduced mass. 
In the adiabatic approximation the total wave function of the 
molecule 



a (r,i?) = cl (fWi?) 



(5) 



is represented as a product of the electronic wave function <p e \ 
and the nuclear wave function \- The electronic wave function 
obeys the clamped nuclei electronic Schrodinger equation 



[H el - E ol (R)] |0 e i) = 0, 



(6) 



while the wave function % is a solution to the nuclear 
Schrodinger equation with the effective potential generated by 
electrons 



V 2 

-T^ + (>eli# n |0el) + E e l(R) - E 3 
A jl n 



|X) = 0, (7) 



where the so called diagonal (or adiabatic) correction 

(0el|-Hn|0cl) = 7, (Vi?^ c l|V_R(/) e l) 
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is a function of R. 

The existence of the electric dipole transitions in HD is due 
to the last term in Eq. (01). This term can be used directly 
as a perturbation. Such an approach is presented in the Ap- 
pendix for a comparison with previous works. In an alterna- 
tive method, inspired by the work of Thorson et al. Il6ll and 
applied here, we introduce a unitary transformation 

H' = U + HU (9) 

to shift the odd term in H a to the potential V in Eq. (0. This 
transformation greatly simplifies further calculations. We 
choose U of the form 

JJ _ e A(ri+f 2 )-V H (1Q) 

with 



Below, we rewrite this matrix element in terms of the elec- 
tric dipole moment function D(R), namely 

D fi = (xf|X>Tl|xi> = <J r f,Mf|fi|Ji,Mi>Ufl, (18) 
D fi = J dRR 2 D(R)x* j! {R)x.h(R), d9) 

D[R) = (UhL^UhL) (20 ) 
1 ' \M B Ma J 

x(4> cX \f-n — v • Vr(V) |0el>. 

The function D(R) depends only on the distance R between 
the nuclei. Although similar, D(R) can not be identified with 
the projection of the dipole moment operator onto the symme- 
try axis, because the direction of R is changed under applied 
unitary transformation. 



2 \M B M a J 



(11) 



and obtain H' with neglecting 0[(m c / /i n ) 2 ] terms, namely 

H' = H + \[H,(? X +r 2 ) ■ Vfl] +0(A 2 ) 

= H - A (n + f 2 ) • V R (V) - — (Vi + Va) • V fi 



//el + 5^ + H' n , 



(12) 



where 



sv -"-r(ik-^) {f ' +f ' } -^ iv) <13) 

V B (V) = i f-Jd + ^£ - % + *°) -* (14, 

2 V r lA r lB r 2A r 2 B J R 



Hi = 



VI (Ea Va)' 



2 /X n 8 jtXn 



(15) 



and n = R/R. 

The El transition between rovibrational levels of the HD 
molecule in the ground electronic state comes now from the 
nonadiabatic correction SV to the electronic potential V. In 
the leading order one uses the adiabatic approximation, and 
the electric dipole moment Dr between some initial 0j and 
final state <f>t is 



D n = (cf> { \f\4>i) 



(16) 



?el — H e l 



■8V\<Wxx) 



1 



Eei — H c \ 



■f\faXl), d7) 



where r = ^2 a r a . We claim, without presenting the proof, 
that the higher order nonadiabatic corrections are smaller by a 
factor of m c / fi n w 10 -3 , and their contribution to can be 
neglected. 



III. NUMERICAL CALCULATIONS 

For the numerical calculation of D(R), the clamped nu- 
clei electronic wave functions were represented in the form of 
linear expansions in the two-electron basis of exponentially 
correlated Gaussian (ECG) functions 



Mrur 2 ) = i(l + A 2 )(l±i) 



(21) 



x exp 



■ ^ Ak.ij(ri — Sfe ) i)(rj s k,j) 



where the matrices and vectors contain nonlinear pa- 
rameters, 5 per basis function, to be variationally optimized. 
The Gaussian centers su were constrained to the internuclear 
axis to preserve the S symmetry. The antisymmetry projector 
(1 + P12) ensures singlet symmetry and the spatial projec- 
tor (1 ± i) — the gerade (+) or ungerade (— ) symmetry with 
respect to inversion in the origin of the coordinate system lo- 
cated at the geometric center of the nuclei. 

The computations were performed independently at 56 in- 
ternuclear distances. In order to check the asymptotic be- 
havior of the dipole moment function, long distances (up to 
R — 12.0 au) were sampled. At every distance R, two 600- 
term basis sets were generated — one, of the symmetry, 
to represent the electronic ground state wave function <f> e \, and 
the other, of the 1 E+ symmetry, to invert the Hamiltonian. 
The parameters of the first basis set were optimized with re- 
spect to the lowest root of the clamped nuclei Hamiltonian 
H c \ and the electronic energy was converged to an accuracy 
of the order of a fraction of microhartree. The nonlinear pa- 
rameters of the second basis were optimized with respect to 
the functional corresponding to the parallel polarizability 



J = 



1 



t>ei\r ■ n 



■ n \(j)e 



(22) 



with the fixed (f> c \ wave function. The basis sets generated 
this way were subsequently employed to evaluate the dipole 



3 



moment D(R), Eq. (l20l i. The proton and the deuteron mass 
used in Eq. ( 1201 were taken from lfl9ll 

M A = M H = 1836.15267247 m e (23) 
M B = Md = 3670.4829654 m e . (24) 

The dipole moment is commonly expressed in the units of de- 
bye (D), and the numerical factor used in this work to convert 
the results from the atomic to debye units was 2.54174623 
D/au. To inspect the saturation of the 1 S+ basis at R = 1.4 
au, we generated an additional 600-term basis set with the 
nonlinear parameters optimized with respect to 

(Mr- VflOO - 1 - v • Vr(V) \<fa) (25) 

and combined this basis set with the original 1 £+ basis. De- 
spite doubling the size of the basis set, the D (R) value has 
changed only on the 8-th significant figure. Hence, we ex- 
pect that all displayed figures of the final result in Table H] are 
significant. 

Numerical values of the D(R) function are presented in Ta- 
ble U and plotted in Figure [TJ For comparison with previous 



TABLE I: D(R) — the electric dipole moment (in 10~ 4 D) as a func- 
tion of the internuclear distance R. All digits are numerically signif- 
icant. 



R/au 


D(R) 


R/a,u 


D(R) 


0.5 


-27.6224 


3.2 


-5.0893 


0.6 


-21.0635 


3.3 


-4.7735 


0.7 


-17.0294 


3.4 


-4.4482 


0.8 


-14.3747 


3.5 


-4.1182 


0.9 


-12.5426 


3.6 


-3.7882 


1.0 


-11.2339 


3.8 


-3.1468 


1.1 


-10.2754 


4.0 


-2.5552 


1.2 


-9.5603 


4.2 


-2.0342 


1.3 


-9.0193 


4.4 


-1.5928 


1.4 


-8.6054 


4.5 


-1.4022 


1.45 


-8.4353 


4.6 


-1.2307 


1.5 


-8.2853 


4.8 


-0.9409 


1.6 


-8.0347 


5.0 


-0.7137 


1.7 


-7.8352 


5.2 


-0.5381 


1.8 


-7.6721 


5.25 


-0.5010 


1.9 


-7.5336 


5.5 


-0.3495 


2.0 


-7.4102 


5.75 


-0.2428 


2.1 


-7.2933 


6.0 


-0.1682 


2.2 


-7.1758 


6.5 


-0.0804 


2.3 


-7.0517 


7.0 


-0.0384 


2.4 


-6.9156 


7.5 


-0.0184 


2.5 


-6.7632 


8.0 


-0.0089 


2.6 


-6.5913 


8.5 


-0.0044 


2.7 


-6.3972 


9.0 


-0.0022 


2.8 


-6.1799 


9.5 


-0.0012 


2.9 


-5.9390 


10.0 


-0.0006 


3.0 


-5.6755 


11.0 


-0.0002 


3.1 


-5.3911 


12.0 


-0.0001 



calculations, the dipole moment function obtained by Ford 
and Browne lfl5ll in the range of R € (0.5, 3.0) au is pre- 
sented in the same Figure Q] The function D(R) behaves as 
R~ 2 at R — ► 0, and as R~ 4 at R —> oo. The singularity at 
R = comes from the neglecting higher order terms in the 
unitary transformation and from the adiabatic approximation. 
At R of the order \J m c /fj, n « 0.03 adiabatic approximation 



- 




R 



FIG. 1: (Color online) Electric dipole moment function D(R) (in 
10 _4 D). Comparison to previous calculations by Ford and Browne 
1 15] (dotted line), n is directed from deuteron to proton. Negative 
D(R) means that electrons are shifted toward the deuteron. 

fails and our formula for D(R) is not valid. At this region, 
however, nuclear wave function is negligible. 

For the calculation of the electric dipole moments, the adi- 
abatic potential of the nuclear Schrodinger equation (0 has 
been composed of the clamped nuclei energy E e \, and the adi- 
abatic correction ((j> c \\H n \(j) c \). For E c \ we used the energy 
points computed with nanohartree accuracy by Cencek from 
1200-term ECG wave functions 112011 . The adiabatic correc- 
tion in Eq. (0 was evaluated by us using Eq. d40l . The adia- 
batic potential curve was then obtained by means of piecewise 
polynomial interpolation. The radial Schrodinger equation 
has been solved numerically using the Le Roy's code 12111 . 
The obtained nuclear wave functions \ of rovibrational levels 
were subsequently used in the evaluation of the dipole transi- 
tion moments of Eq. (fl~8l for the J — > J + 1 (branch R) and 
for the J — > J — 1 transitions (branch P). 



IV. RESULTS AND DISCUSSION 



Our electric dipole moments for the transitions between the 
lowest vibrational and rotational levels are listed in Table [TT] 
Except for 0-0 transitions, they are in good agreement with the 
previous calculation by Thorson et al. For the lowest band all 
theoretical predictions differ slightly from each other. Our re- 
sults are numerically accurate to all digits shown, but the last 
digit is uncertain due to the neglected 0(m c //x n ) pa 8 x 10~ 4 
higher order nonadiabatic corrections. These corrections have 
been also neglected in calculations of lfl5l [l6Tl . so in prin- 
ciple these calculations should agree with each other. Con- 
sidering calculations presented in [14], we note that in the 
initial expression for D^, Wolniewicz uses H'^, Eq. d26t . as 
a perturbation, and assumes the adiabatic approximation for 
the wave function, but in the denominator he includes H' n 
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TABLE II: Experimental and theoretical electric dipole transition moments Da (in 10 4 D). Relative uncertainty of our results due to the 
nonadiabatic corrections to Dti is about 10 -3 . 



Reference 


P(3) 


P(2) 


P(l) 


R(0) 


R(l) 


R{2) 


R(3) 


experiment [3] 






0-0 




9.36(30) 


8.00(20) 


9.79(30) 


experiment [4] 










8.78(2) 


8.47(2) 


10.21(2) 


theory [44] 








8.36 


8.38 


8.39 


8.41 


theory [45J 


8.282 


8.297 


8.306 


8.306 


8.297 


8.282 


8.262 


theory 1161 


8.440 


8.455 


8.463 


8.463 


8.455 


8.440 


8.420 


this work 


8.536 


8.551 


8.560 


8.560 


8.551 


8.536 


8.516 



1-0 



experiment [6] 


0.330(40) 


0.405(30) 


0.450(30) 


0.515(20) 


0.550(30) 


0.615(30) 


0.655(40) 


experiment [7] 


0.340(22) 


0.379(12) 


0.435(11) 


0.504(12) 


0.533(14) 




0.609(13) 


theory [44] 








0.598 


0.628 


0.656 


0.685 


theory [45] 


0.401 


0.445 


0.485 


0.560 


0.594 


0.623 


0.650 


theory 1161 


0.374 


0.421 


0.466 


0.552 


0.592 


0.630 


0.665 


this work 


0.3776 


0.4248 


0.4708 


0.5579 


0.5983 


0.6362 


0.6714 



2-0 



experiment [6] 






0.17(2) 


0.19(2) 


0.20(2) 






theory [44] 








0.160 


0.166 


0.170 


0.174 


theory [45J 


0.156 


0.167 


0.176 


0.192 


0.199 


0.206 


0.210 


theory 1161 


0.156 


0.167 


0.179 


0.200 


0.210 


0.219 


0.228 


this work 


0.1576 


0.1692 


0.1805 


0.2022 


0.2122 


0.2216 


0.2301 



3-0 



experiment [61 
theory [44] 
theory [15J 
theory 1161 
this work 



0.068 

0.0698 

0.0705 



0.072 

0.0742 

0.0749 



0.076 

0.0786 

0.0794 



0.0795(35) 

0.100 

0.082 

0.0870 

0.0878 



0.0800(50) 

0.102 

0.084 

0.0909 

0.0918 



0.103 
0.085 
0.0945 
0.0955 



0.104 
0.087 
0.0979 
0.0989 



experiment [8] 
theory [44] 
theory [45J 
theory 1161 
this work 



0.033 

0.0324 

0.0327 



0.035 

0.0345 

0.0348 



4-0 

0.0397(26) 

0.038 

0.0365 

0.0369 



0.0417(24) 


0.0425(21) 


0.0459(26) 


0.0514(53) 


0.056 


0.056 


0.055 


0.053 


0.039 


0.040 


0.042 


0.042 


0.0405 


0.0425 


0.0442 


0.0458 


0.0409 


0.0428 


0.0446 


0.0462 



5-0 



experiment [91 


0.0105(25) 


0.0124(21) 


0.0143(17) 


0.0181(17) 


0.0200(21) 


0.0219(25) 





experiment [81 








0.0207(20) 


0.0214(14) 


0.0231(21) 




theory [15J 


0.020 


0.021 


0.021 


0.023 


0.023 


0.024 


0. 


theory 1161 


0.0163 


0.0173 


0.0184 


0.0205 


0.0215 


0.0225 





this work 


0.0164 


0.0175 


0.0186 


0.0207 


0.0217 


0.0227 






0.0238(29) 



from Eq. dl5l l. This expression, in comparison to our, in- 
volves some higher order terms, namely X± of Eq. (l37l i from 
Appendix. However, we show in the Appendix the cancella- 
tion of significant contributions involving the second deriva- 
tive of x between X\ and the other higher order contributions 
Xi (Eq. (l38l >) and X% (Eq. (|39i>), which has been neglected 
in Wolniewicz calculations by assuming the adiabatic wave 
function. Therefore, we think, slight difference with results of 
Wolniewicz in ffl may come from less consistent treatment 
of higher order nonadiabatic effects. 

In comparison to experimental values, we observe a mod- 
erate agreement for all transitions, but the 0-0 ones. Here, 



our results, as well as the other theoretical calculations, dif- 
fer from the experiment by several standard deviations. We 
note however, that the measurements are most cumbersome 
for these transitions. As a consequence, experimental values 
significantly change with the rotational number J for v = 0, 
which can not be justified by theoretical analysis. Within the 
ground vibronic state, the nuclear wave functions correspond- 
ing to the lowest rotational levels are localized near the aver- 
age internuclear distance R , and differ very little from each 
other. For this reason subsequent transition moments must 
change slowly with the rotational quantum number J, and are 
approximately equal to D(R Q ), but the experimental results 
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of JH01 are not consistent with theoretical predictions. 

V. CONCLUSIONS 

We have derived a simple expression for the electric dipole 
transition rates, in terms of the dipole moment function D(R), 
and performed precise calculations of D(R) in a wide range 
of R. The obtained formula, can easily be extended to other 
diatomic molecules, consisting of two isotopes of the same 
element. Our results for the dipole moments of HD molecule 
are numerically accurate to four digits, in moderate agreement 



with previous calculations in lfl4l[T5i[T6ll and experimental re- 
sults of EH SUBS Hi ( see TableHD. We estimate, that the 
relativistic and nonadiabatic corrections are of relative order 
of 10~ 4 and 10~ 3 , correspondingly. As no other effect may 
alter the theoretical predictions, we suppose, that our results 
are even more accurate than the experimental values obtained 
so far. 
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Appendix 

For the comparison with previous works Ihil [l5tl . which 
used the term 



(26) 



as a perturbation, we derive the formula for Dr, Eq. (ill , 
using the nonadiabatic perturbation theory. For this, one has 
to abandon the assumption in Eq. (0 of a separation of the 
electronic wave function from the nuclear one. Namely, the 
total wave function 



4>e\ X + <50n 



(27) 



will be the sum of the adiabatic solution and a nonadiabatic 
correction. The nonadiabatic correction <50 na is decomposed 
into two parts 



<5</>na = 4>cl 5x + 8'(f>na,, 

which obey the following orthogonality conditions 

(<5'0na|</>el)cl = 0, 

(Sx\x) = 0, 



(28) 



(29) 
(30) 



with the normalization \4>\4> a ) = 1- In the leading order 
of perturbative treatment, the nonadiabatic corrections to the 
wave function are the following [18] 

\8'<t>$) = , F * ffnl^dX), (3D 



\Sx) = 



[J3„ — E e \ — H n — (-Hn/eJ 



-((f> c \\H n \S'(j) na ) cV 



(32) 



where the prime in the denominator denotes subtraction of the 
reference state from the Hamiltonian inversion. For the calcu- 
lation of D(R) one needs also the second order correction 



\5'€1) = 



x [H n + E cl - E a 



1 



{E c l — Hel) 



7 H n \4> e ix), 



(33) 



where 8% is given by Eq. 

The derivation of the formula <f2Qb proceeds as follows. One 
takes Eq. ( TToT i with perturbed wave functions, 



-Dfi = (0el (Xf + <5Xf) + <5'</>f,naH</>cl (Xi + <>Xi) + <*'0i,na), 

(34) 



4 1} + 4 2) 



i 



#fi' = W>ei»|r 

£/ H — rl, 



belXf\K 



I — flel 
1 



E e l — H c i 



r\^>ei Xi) , 



(35) 
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and is given in Eq. (|44| |. In the one separates out 
the electronic matrix elements from the nuclear ones, namely 



given by 



-2{M- B -M- A )J d ' R {** 
x2 (^ci|r fc 1 V fi ■ V V a |0 e i) 

a 

+(xjf v£xi + xi vixf) 

X (^lk fc -^^- E^l^l)}. (36) 



where the superscripts j and k are the Cartesian indices. The 

(l)k 

second term in braces is integrated by parts and D\ ' be- 
comes 



D 



(2) 



PelXi\r 
1 



(E e \ — H e \) 



7 [H n + E c \ — Ei ta \ 



(E c \ - H c \) 

+{ ^ H »(E cl -H cl y f (E cl -H cl y Hn 

x|0elXi) + {faXiWn 



(E el - H el y 



; [Hn + E e l — £f,a] 



1 



(Eel - H c 



y ^l^el Xi)- 



(44) 



D 



(i)k _ 



where 



x[(^|V^)-<V^f|^>], 



where we neglected <$x- For low lying rovibrational states 
(37) $x/x is small, namely of 0(m a /fi n ) and thus its magnitude 
is of order 10~ 3 , if not less, and thus negligible. H n is de- 
composed into the even H' n , Eq. ( TTBT l. and the odd H[[ parts 
Eq. d26l ). involves the terms with a single power of H' n . 

(38) 

The resulting 6 terms in Z)g 2 ^ we group into pairs as follows 



= -(H c \- Eeijrrier 1 \<f) c \). 

One takes of the Schrodinger equation (|6]l to obtain 
1 



(E C 1 - Hel) 
1 

Ed — H c \ 

+r k 
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(39) 



(40) 



1 



Ee\ — H c l 



(Eel — H e l) 



(41) 



Vr<& = -V 3 R (V- Earner* fa 



+(E e i - H e i)m c r 3 



(Eel — H e l) 



(42) 



The gradient of the electronic functions with respect to the 
internuclear distance in Eq. (f37T > is replaced by Eqs. (HTT i and 
(l42l . Among the four terms, two cancel out and the two other 
are the same, so the transition dipole moment takes the form 



= ^1+^2+^3, 



(45) 



Xl = (0elXf|r 



1 



(Eel — H e 



[H' n + Ed — Ei ta ] 



x ... 1 , HZ\faXi> + {<Uxt\K 1 



(Ed — H c 
[H n + Ed — Ef >a ] 



(E e l — Hgl) 



(E e l — H e 



■r-^eiXi), (46) 



X 2 = ((t>eiXf\H' nT = — ^—=-^r 1 H' n ' \(f>elXi) 
+ (<l)oiXi\H' n ' ,„ 1 „ : r— — ~ H' n \<f> e iXi), 

(±j e l — Hcl) [Eel — Hel) 

(47) 

X 3 = (</>eiXf|r-7- TT~\ H ^T^ 5-17^1^1X1) 

(£j e \ — £t e l) (tiel — -tlcl) 

+ (<t>eiXf\H n T = 1 - E -^ K — 1 - r r|0 el Xi). 



(E c l — H c i)' 11 (Ed ~ H e l) 



(48) 



Ji)k _ ( me_ _ ro e \ 



D ="~ = ug-sjy^* (43) 

x(0ei|r fe - 1 - r ■ V R (V - Eei) [fa). 

hi e l — Hcl 

One notes that it differs from Eq. < f20b only by the presence of 

~~* (2)k 

VrE c \. We show below, that this term cancels out with D & 



H' n involves two derivatives over R and H[[ a single deriva- 
tive. Consider terms with three derivatives of x- Since x sat- 
isfies Eq. ©, the second derivative of x coming from H' n can- 
cels with E a — E c i(R), leaving a small term (H n ) c i and the 
derivative of E c i- Thus no term with three derivatives of x 
is present. We will show below that no terms involving any 
derivatives of x are present. Each X; includes two derivatives 
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terms 

+* (V' B Vi xf )] (&,|r* ' r'lV^ei), 

(49) 

+(Vixf) (V,, X,)] tfal^r* - 1 |Vi,fe), 

iicl — _H i 

(50) 

X* « -X*-X 2 fc , (51) 



but all of them cancel out in the sum. A single derivative of 
X has to be of the form Vn(Xf Xi)> what can be integrated by 
parts and transformed into derivatives of c i and the resolvent. 
Considering terms without derivatives of \ all of them are of 
order 0[(m e /m n ) 2 ] but one, which involves VjiE e \ which 
arised from the commutator [H' n + E C \(R) — E a , V«] in X\ 
in Eq. (46). It is of the form 

5 «' " *«*-^ 5 7r.V H ( £ .,)| tel > 

(52) 

which together with Dg gives the leading correction to the 
transition dipole moment Eq. ( fT71 ). 



